UMN-TH-2617/07 
FTPI-MINN-07/27 
September 2007 



Colliders and Cosmology 



o 
o 

(N 

Or 
<D 

o 

a,; 

in: 



> 

m 
o 
m 
cn 

ON 

o 
o 



X 



Keith A. Olive 1 a 

William I. Fine Theoretical Physics Institute, University of Minnesota, Minneapolis, MN 55455, USA 



Abstract. Dark matter in variations of constrained minimal supersymmetric standard models will 
be discussed. Particular attention will be given to the comparison between accelerator and direct 
detection constraints. 
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1 Introduction 



Evidence for dark matter in the universe is available 
from a wide range of observational data. In addition 
to the classic evidence from galactic rotation curves 
[1], the analysis of the cosmic microwave background 
anisotropics leads to the conclusion that the curvature 
of the universe is close to zero indicating that the sum 
of the fractions of critical density, £2, in matter and a 
cosmological constant (or dark energy) is close to one 
[2]. When combined with a variety of data including 
results from the analysis of type la supernovae obser- 
vations [3] and baryon acoustic oscillations [I] one is 
led to the concordance model where Q m ~ 0.23 and 
Qa ~ 0.73 with the remainder (leading to fttot = 1) 
in baryonic matter. Other dramatic pieces of evidence 
can be found in combinations of X-ray observations 
and weak lensing showing the superposition of dark 
matter (from lensing) and ordinary matter from X-ray 
gas [5] and from the separation of baryonic and dark 
matter after the collision of two galaxies as seen in the 
Bullet cluster [BJ. For a more complete discussion see 

Here, I will adopt the results of the three-year data 
from WMAP [2] which has determined many cosmo- 
logical parameters to unprecedented precision. Of par- 
ticular interest is the determination of the total matter 
density (relative to the critical density), f2 m h 2 and the 
baryonic density, fl\>b? . In the context of the vlCDM 
model, the WMAP only results indicate 
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The difference corresponds to the requisite dark mat- 
ter density 
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or a 2cr range of 0.0882 - 0.1204 for f2cDAih 2 
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2 Constrained MSSM models 

In its generality, the minimal supersymmetric stan- 
dard model (MSSM) has over 100 undetermined pa- 
rameters. But in addition to relieving the helplessness 
of an analysis with so many free parameters, there are 
good arguments based on grand unification [5] and su- 
pergravity [5] which lead to a strong reduction in the 
number of parameters. I will assume several unification 
conditions placed on the supersymmetric parameters. 
In all models considered, the gaugino masses are as- 
sumed to be unified at the GUT scale with value, m^ 2 , 
as are the trilinear couplings with value Aq. Also com- 
mon to all models considered here is the unification 
of all soft scalar masses set equal to mo at the GUT 
scale. With this set of boundary conditions at the GUT 
scale, we can use the the radiative electroweak symme- 
try breaking conditions by specifying the ratio of the 
two Higgs vacuum expectation values, tan/3, and the 
mass, Mz, to predict the values of the Higgs mixing 
mass parameter, fi and the bilinear coupling, B. The 
sign of {j, remains free. This class of models is often 
referred to as the constrained MSSM (CMSSM) [H 

Although the CMSSM is often confused with in 
name with mSUGRA, i.e. models based on minimal su- 
pergravity, [15] , the latter employ two additional con- 
straints [15] . In addition to the conditions listed above, 
these are: a relation between the bilinear and trilinear 
couplings Bq = Aq — mo, and the relation between 
the gravitino mass and soft scalar masses, m 3 / 2 = mo. 
When electroweak symmetry breaking boundary con- 
ditions are applied, this theory contains only mm, mo, 
and Aq in addition to the sign of fj,, as free parameters. 
The magnitude of fi as well as tan /3 are predicted. 

I will also discuss a less constrained model, the 
NUHM, in which the Higgs soft masses are not uni- 
fied at the GUT scale [T71I18] . In this class of models, 
both fi and the Higgs pseudo scalar mass become free 
parameters. 

In all of the models discussed below, I will assume 
unbroken R — parity and and hence the lightest super- 
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symmetric particle is stable. This will often, but not 
always be the neutralino [19] . 



3 The CMSSM after WMAP 

For a given value of tan/3, A , and sgn(fi), the result- 
ing regions of acceptable relic density and which satisfy 
the phenomenological constraints can be displayed on 
the mi/2 — m plane. In Fig. [IJl, the light shaded region 
corresponds to that portion of the CMSSM plane with 
tan /3 = 10, A = 0, and /i > such that the computed 
relic density yields the WMAP value given in eq. ([2]) 
[13] . The bulk region at relatively low values of 
and Too, tapers off as my 2 is increased. At higher val- 
ues of mo, annihilation cross sections are too small to 
maintain an acceptable relic density and f2 x h 2 is too 
large. Although sfermion masses are also enhanced at 
large mi/2 (due to RGE running), co-annihilation pro- 
cesses between the LSP and the next lightest sparti- 
cle (in this case the f ) enhance the annihilation cross 
section and reduce the relic density. This occurs when 
the LSP and NLSP are nearly degenerate in mass. The 
dark shaded region has < m x and is excluded. The 
effect of coannihilations is to create an allowed band 
about 25-50 GeV wide in mo for my 2 950 GeV, or 
m i/2 400 GeV, which tracks above the — m x 
contour [20] . 

Also shown in Fig. [T^ are the relevant phenomeno- 
logical constraints. These include the LEP limits on 
the chargino mass: m x ± > 104 GeV [21] , on the selec- 
tron mass: mg > 99 GeV [22] and on the Higgs mass: 
m/j > 114 GeV [23]. The former two constrain mi/2 
and mo directly via the sparticle masses, and the lat- 
ter indirectly via the sensitivity of radiative corrections 
to the Higgs mass to the sparticle masses, principally 
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FeynHiggs [23] is used for the calculation of rrih- 
The Higgs limit imposes important constraints princi- 
pally on mi/2 particularly at low tan/3. Another con- 
straint is the requirement that the branching ratio for 
b — ► 57 is consistent with the experimental measure- 
ments [25j . These measurements agree with the Stan- 
dard Model, and therefore provide bounds on MSSM 
particles [26], such as the chargino and charged Higgs 
masses, in particular. Typically, the b — > sj constraint 
is more important for /i < 0, but it is also relevant 
for /i > 0, particularly when tan/3 is large. The con- 
straint imposed by measurements of b — > sj also ex- 
cludes small values of my 2 . Finally, there are regions 
of the (mi/2,mo) plane that are favoured by the BNL 
measurement [27] of g^ — 2 at the 2-cr level, correspond- 
ing to a deviation from the Standard Model calcula- 
tion [28] using e + e" data. 

Another mechanism for extending the allowed re- 
gions in the CMSSM to large m x is rapid annihilation 
via a direct-channel pole when m x ~ \rriA [TUlll2| . 
Since the heavy scalar and pseudoscalar Higgs masses 
decrease as tan/3 increases, eventually 2m x ~ tua 
yielding a 'funnel' extending to large m^ 2 and mo at 
large tan /3, as seen in Fig. [T)d. As one can see, the im- 
pact of the Higgs mass constraint is reduced (relative 
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Fig. 1. The (w.i/2,mo) planes for (a) tan/3 = 10 and 

> 0, assuming Ao = 0, m t = 175 GeV and mb{mb)'sM ~ 
4.25 GeV. The near-vertical (red) dot-dashed lines are the 
contours mi, = 114 GeV, and the near-vertical (black) 
dashed line is the contour m x ± = 104 GeV. Also shown 
by the dot-dashed curve in the lower left is the corner ex- 
cluded by the LEP bound of mg > 99 Ge V. The medium 
(dark green) shaded region is excluded by b — ^ s'y, and the 
light (turquoise) shaded area is the cosmologically preferred 
region. In the dark (brick red) shaded region, the LSP is the 
charged fi . The region allowed by the E821 measurement 
of at the 2-o level, ts shaded (pink) and bounded by solid 
black lines, with dashed lines indicating the 1-a ranges. In 
(b), tan/3 = 50. 



to the case with tan/3 = 10) while that of b — > sj is 
enhanced. 

Shown in Fig. [2] are the WMAP lines [13 of the 
( TO i/2i m o) plane for \x > and values of tan/3 from 
5 to 55, in steps zi(tan/3) = 5. We notice immedi- 
ately that the strips are considerably narrower than 
the spacing between them, though any intermediate 
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Fig. 2. The strips display the regions of the ^^^o) 
plane that are compatible with the WMAP determination of 
fl x h 2 and the laboratory constraints for fi > and tan/3 = 
5, 10, 15, 20, 25, 30, 35, 40, 45, 50, 55. The parts of the strips 
compatible with g M — 2 at the 2-a level have darker shading. 



point in the (m^ 2 ,mo) plane would be compatible 
with some intermediate value of tan/3. The right (left) 
ends of the strips correspond to the maximal (mini- 
mal) allowed values of mi/2 and hence m x . The lower 
bounds on 7774/2 are due to the Higgs mass constraint 
for tan/3 < 23, but are determined by the b — > S7 
constraint for higher values of tan (3. 

Finally, there is one additional region of accept- 
able relic density known as the focus-point region [29 , 
which is found at very high values of uiq. An exam- 
ple showing this region is found in Fig. [31 plotted for 
tan/3 = 10, /j > 0, and m t — 175 TeV. As m is 
increased, the solution for \i at low energies as deter- 
mined by the electroweak symmetry breaking condi- 
tions eventually begins to drop. When /j < m 1 / 2 , the 
composition of the LSP gains a strong Higgsino com- 
ponent and as such the relic density begins to drop 
precipitously. As mo is increased further, there are no 
longer any solutions for fi. This occurs in the shaded 
region in the upper left corner of Fig. [3] 



4 Direct detection 

Direct detection techniques rely on the neutralino nu- 
cleon scattering cross-section. In general, there are two 
contributions to the low-energy effective four-fermion 
Lagrangian which are not velocity dependent. These 
can be expressed as spin-dependent and scalar or spin- 
independent interactions. Here I will discuss only spin- 
independent elastic x-nucleon scattering stemming from 



Fig. 3. As in Fig. QJs, where the range in mo is extended 
to 5 TeV. In the shaded region at very high mo, there are 
no solutions for fi which respect the low energy electroweak 
symmetry breaking conditions. 



type quarks (i — 2). The expression for a^i can be 
found in [30] . 

The scalar cross section can be written in terms of 

<j 3 = M [ Z f p + (A — Z)f n } 2 , 



C = 



(3) 



which is to be summed over the quark flavours q, and 
the subscript i labels up- type quarks (i — 1) and down- 



(4) 



where m r is the reduced LSP mass, 
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the parameters ff^ are defined by 

m pfTq = (p\ m qQl\p) = m g B q , 



(0) 



fro = 1 ^ T, q =u,d,s frq ED, and /„ has a similar ex- 
pression. The needed matrix elements are determined 
in part by the 7T-nucleon £ term, which is given by 

ovjv = £ = ^(m u + m d )(B u + B d ). (7) 

and the strangeness contribution to the proton mass, 
2B S 



Ap) 



y 



B„ + Ba 



(8) 



where 00 is the change in the nucleon mass due to the 
non-zero u, d quark masses, which is estimated on the 
basis of octet baryon mass differences to be tro = 36 ±7 
MeV 32 j. Values of the matrix elements are given in 
the table for 3 choices of £ and hence y. As one can 
see, we expect there will be a strong dependence of the 
cross section on the proton strangeness, y. 
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Table 1. Matrix elements in terms of proton strangeness the LEP m n constraint for tan/3 = 10 and the b 

constraint for tan/3 = 50. 



S7 



E (MeV) 


y 






St. 


45 


0.2 


0.020 


0.026 


0.117 


64 


0.44 


0.027 


0.039 


0.363 


36 





0.016 


0.020 


0.0 



5 Indirect sensitivities 

Measurements of electroweak precision observables 
(EWPO) as well as B-physics observables (BPO) can 
provide interesting indirect information about the su- 
persymmetric parameter space. We have already seen 
the impact of measurements of the anomolous mag- 
netic moment of the muon, the branching raio of b — > 
S7 in addition to the non-discovery of charginos and 
the Higgs boson at LEP which impose significant lower 
bounds on m^- 

In [33 l l34 j l35j . we considered the following EWPO: 
the W boson mass, Mw, the effective weak mixing an- 
gle at the Z boson resonance, sin 2 9 e g, the width of the 
Z, Fz, the Higgs mass, rrih, and (g^ — 2) in addtion 
to the BPO: b decays b — > 37, B s — > B u — > tv t 

and the J5 S mass mixing parameter AMb b - We per- 
formed the analysis of the sensitivity to TO1/2 moving 
along the WMAP strips with fixed values of A§ and 
tan/3. The experimental central values, the present ex- 
perimental errors and theoretical uncertainties are as 
described in [35] . Assuming that the nine observables 
listed above are uncorrelated, 
formed with 



a x nt has been per- 
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Here R n xp denotes the experimental central value of 
the nth observable (Mw, sin 2 9 e s, Z~!z, {g — 2) M and 
BR(6 -> sj), BR(B U -> rz/ T ), AM Bs ), R n hco is the cor- 
responding MSSM prediction and a n denotes the com- 
bined error. Additionally, cr™ ln is the minimum com- 
bined error over the parameter space of each data set, 
and Xm h an d Xb denote the \ 2 contribution coming 
from the experimental limits on the lightest MSSM 
Higgs boson mass and on BR(£> S — > respec- 
tively [33] , 

Details of the analysis [35] were discussed in [36] 
[37] and here I show the net result of the \ 2 analysis 
in the CMSSM in Fig. |4] using the combined \ 2 values 
for the EWPO and BPO, computed from Eq. M We 
see that the global minimum of \ 2 ~ 4.5 for both val- 
ues of tan p. This is quite a good fit for the number 
of experimental observables being fitted. For both val- 
ues of tan/3, the focus-point region is disfavoured by 
comparison with the coannihilation region, though this 
effect is less important for tan/3 = 50. For tan/3 = 10, 
m i/2 ~ 300 GeV and Aq > are preferred, whereas, 
for tan/3 = 50, m 1 / 2 ~ 600 GeV and A < are pre- 
ferred. This change-over is largely due to the impact of 



14 



12 



10 



CO 
"D 
O 




• • • 

y 

. CMSSM, ir > 0, m, = 171.4 GeV - 
,•* . tan(i = 10, A = 

• tan|3 = 10, A = +m 1/2 
o tanp = 10, A = -m 1/2 
. tanp = 10, A = +2 m 1;2 
tanp = 10, A = -2 m 




1000 



2 - 



lJ_ 



CMSSM, n > 0, m t = 1 71 .4 GeV- 
tanp = 50, A = 
tanp = 50, A = +m 1(2 
tanp = 50, A = -m 1;2 
tanp = 50, A = +2 m y2 
tanp = 50, A n = -2 m„„ 



200 400 600 800 1000 1200 1400 1600 1800 2000 

m 1/2 [GeV] 

Fig. 4. The combined \ 2 function for the electroweak ob- 
servables Mw, sin 2 9 e s, Fz, (g — 2) M , rrih, and the b physics 
observables BR(6 -> sj), BR(B S -> BR(B u -» 

tv t ) and AM Bs , evaluated in the CMSSM for tan/3 = 10 
(a) and tan/3 = 50 (b) for various discrete values of Ao- 
We use mt = 171.4 ± 2.1 GeV and m b {m b ) = 4.25 ± 0.11 
GeV, and mo is chosen to yield the central value of the 
cold dark matter density indicated by WMAP and other 
observations for the central values of mt and mj(mi). 



6 Direct Detection in the CMSSM 

In Fig. 03 we display the expected ranges of the spin- 
independent cross sections in the CMSSM when we 
sample randomly tan/3 as well as the other CMSSM 
parameters |38] . All points shown satisfy all phenomeno- 
logical constraints and have fi x h 2 less than the WMAP 
upper limit. Since models with low S7 x h 2 can not be 
excluded if another form of dark matter is assumed to 
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exist, these models have been included, but their cross 
sections have been scaled by a factor Q X J Ocdm- in 
Fig. [5^, E-xN — 45 MeV has been chosen, and in [SJd 
results for S^n = 64 MeV are shown for comparison. 
As one can see, the cross sections shift higher by a 
factor of a few for the larger value of H^n- 

Also shown on the plot are the current CDMS [35] 
and XENON10 [10] exclusion curves which place an 
upper limit on the scattering cross section. As one can 
see, the current limits have only just now begun to 
probe CMSSM models. 
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Fig. 5. Scatter plots of the spin-independent elastic- 
scattering cross section predicted in the CMSSM for (a) 
E nN = 45 MeV and (b) 64 MeV. 



7 GUT-less models 

The input scale at which universality is assumed in 
CMSSM models is usually taken to be the SUSY GUT 



scale, Mqut ~ 2 x 10 16 GeV. However, it may be more 
appropriate in some models to assume the soft SUSY- 
breaking parameters to be universal at some different 
input scale, Mj„. Specific scenarios in which the soft 
SUSY-breaking parameters may be universal at a scale 
below Mqut occur in models with mixed modulus- 
anomaly mediated SUSY breaking, also called mirage- 
mediation [JT] , and models with warped extra dimen- 
sions 42J. In the case of mirage- mediation, the uni- 
versality scale is the mirage messenger scale, which is 
predicted to be M m ~ 10 10 - 10 12 GeV in the case of 
KKLT-style moduli stabilization [43] . In other models, 
the universality scale may lie anywhere between 1 TeV 
and Mpi. 

In the CMSSM with universality imposed at the 
GUT scale, the one-loop renormalizations of the gaug- 
ino masses M a , where a = 1,2,3, are the same as 
those for the corresponding gauge couplings, a a . Thus, 
at the one-loop level the gaugino masses at any scale 
Q < Mgut can be expressed as 

M a (Q) = ""j Q) M Mgot), (10) 

where M a (MGUT) = m \/2- O n the other hand, in a 
GUT-less CMSSM, where the gauge-coupling strengths 
run at all scales below the GUT scale but the soft 
SUSY-breaking parameters run only below the lower 
universality scale, Mi n , at which all the gaugino masses 
are assumed to be equal to TO1/2 = M a (Mi n ), we have 

at the one-loop level. Since the runnings of the cou- 
pling strengths in GUT and GUT-less CMSSM scenar- 
ios are identical, the low-energy effective soft gaugino 
masses, M a (Q), in GUT-less cases are less separated 
and closer to mx/2 than in the usual GUT CMSSM. 
Similarly, the sfermion masses also run less and have 
masses closer to mo at the weak scale than their cor- 
responding masses in the CMSSM. 

One of the most dramatic changes when Mi n is 
lowered from the GUT scale is its effect on the foot- 
print in the (wi/2, Wo) plane of the constraint on the 
relic abundance of neutralinos inferred from WMAP. 
In the GUT-less CMSSM scenario 0J, as the univer- 
sality scale is lowered to M\ n ~ 10 12 GeV, the co- 
annihilation strip, the funnel and the focus point re- 
gions of the CMSSM approach each other and merge, 
forming a small WMAP-preferred island in a sea of 
parameter space where the neutralino relic density is 
too small to provide all the cold dark matter wanted 
by WMAP. 

As found in [44], and discussed in [45], there are 
already changes as the universality scale is lowered 
to M in = 10 14 GeV. The allowed focus-point region 
starts to separate from the LEP chargino bound, mov- 
ing to larger m^ 2 . For Mj„ = 10 13 GeV, we found that 
the allowed focus-point region also dips further down, 
away from the electroweak vacuum condition bound- 
ary, while the coannihilation strip moves up and far- 
ther away from the region where the stau is the LSP. 
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Fig. 6. Examples of (m 1( / 2 ,mo) planes with tan/3 = 10 
and Aq = but with different values of Mi n . (a) Mi n = 
10 12 ' 5 GeV, (b) M in = 10 12 GeV. In each panel, we show 
contours representing the LEP lower limits on the chargino 
mass (black dashed line), a Higgs mass of 114 GeV (red 
dashed), and the more exact (and relaxed) Higgs bound (red 
dot-dashed). We also show the region ruled out because the 
LSP would be charged (dark red shading), and that excluded 
by the electroweak vacuum condition (dark pink shading). 
The region favoured by WMAP has light turquoise shading, 
and the region suggested by — 2 at 2-a has medium (pink) 
shading, with the 1-a contours shown as black dashed lines. 



shown in panel (a). We now see clearly two distinct 
regions of the plane excluded due to an excess relic 
density of neutralinos; the area enclosed by the atoll 
and the slice between the lower funnel wall and the 
boundary of the already-excluded r-LSP region. 

In panel (b) of of Fig. © for M in = 10 12 GeV, 
the focus-point and coannihilation regions are fully 
combined and the atoll has mostly filled in to be- 
come a small island of acceptable relic density. To the 
right of this island is a strip that is provided by the 
lower funnel wall. The strip curves slightly as m\n in- 
creases then takes a sharp plunge back down towards 
the boundary of the region where the stau is the LSP, 
a feature associated with the \X h + A threshold. 
Reduction in the universality scale from this point re- 
sults in the lower funnel wall being pushed down into 
the excluded r LSP region and total evaporation of 
the island as seen in Fig. [7] and only a small residual 
turquoise region at large m\i% where the relic density 
is within the WMAP limits. At all other points in the 
visible part of the (mi/2,mo) plane the relic density 
of neutralinos is too low to provide fully the cold dark 
matter density preferred by WMAP et al. Of course, 
these SUSY models would not be excluded if there is 
another source of cold dark matter in the universe. 
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Fig. 7. As in Fig. [H/or M in = 10 10 GeV. 



Another remarkable feature at this value of M; ra is the 
appearance of the rapid-annihilation funnel, familiar 
in the GUT-scale CMSSM at large tan/3, but an un- 
familiar feature for tan/3 =10. 

As the universality scale is further decreased to 
M in = 10 12 - 5 GeV, as shown in panel (a) of Fig. [TJ1 
the atoll formed by the conjunction of what had been 
the focus-point and coannihilation strips has shrunk, 
so that it lies entirely within the range of {m 1 / 2 ,mo) 



8 mSUGRA models 

As discussed earlier, mSUGRA models, in contrast to 
the CMSSM, have an additional boundary condition, 
namely the value of the bilinear £?-term is fixed at 
the GUT scale to B = A — niQ. As a consequence, 
one is no longer free to choose a fixed value of tan/3 
across the m\/2 — m o plane. Instead, tan/3 is derived 
for each value of rai/2)TOo an d [1 6j . Phenomeno- 
logically distinct planes are rather determined by a 
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choice for Ao/my. In Fig.0 two such planes are shown 
with (a) Ao/mo = 3 — V3 as predicted in the simplest 
model of supersymmetry breaking [55] and with (b) 
Ao/ma — 2.0. Shown are the contours of tan/3 (solid 
blue lines) in the (m 1 / 2 , mo) planes. Also shown are the 
contours where m x ± > 104 GeV (near-vertical black 
dashed lines) and rrih > 114 GeV (diagonal red dash- 
dotted lines). The regions excluded by b — ► 57 have 
medium (green) shading, and those where the relic 
density of neutralinos lies within the WMAP range 
have light (turquoise) shading. The region suggested 
by ~ 2 at 2-cr has very light (pink) shading. As 
one can see, relatively low values of tan (3 are obtained 
across the planes. 

Another difference between the CMSSM and mod- 
els based on mSUGRA concerns the mass of the grav- 
itino. In the CMSSM, it is not specified and and can 
be taken suitably large so that the neutralino is the 
LSP (outside the stau LSP region). In mSUGRA, the 
scalar masses, mo are determined by (and equal to) the 
gravitino mass. In Fig. [5J the gravitino LSP and the 
neutralino LSP regions are separated by dark (choco- 
late) solid lines. Above this line, the neutralino (or 
stau) is the LSP, whilst below it, the gravitino is the 
LSP 47J. As one can see by comparing the two panels, 
the potential for neutralino dark matter in mSUGRA 
models is dependent on Ag/mo- In panel (a), the only 
areas where the neutralino density are not too large oc- 
cur where the Higgs mass is far too small or, at higher 
mo the chargino mass is too small. At larger Ag/nia, 
the co-annihilation strip rises above the neutralino- 
gravitino delineation. In panel (b), we see the familar 
co-annihilation strip. It should be noted that the focus 
point region is not realized in mSUGRA models as the 
value of /i does not decrease with increasing mo when 
Ao/rriQ is fixed and Bq = Aq — mo- There are also no 
funnel regions, as tan/3 is never sufficiently high. 

In the gravitino LSP regions, the next-to-lightest 
supersymmetric particle (NSP) may be either the neu- 
tralino or stau which is now unstable. The relic den- 
sity of gravitinos is acceptably low only below the 
dashed (pink) line. This excludes a supplementary do- 
main of the (jtoi/2, Too) plane in panel (a). However, 
the strongest constraint is provided by the effect of 
neutralino or stau decays on big bang neucleosynthe- 
sis (BBN) [58]. Outside the light (yellow) shaded re- 
gion, the decays spoil the success of BBN. Note that 
in panel (b), there remains a region which is excluded 
because the stau is the LSP. 

Recently, new attention has been focussed on the 
regions with a stau NSP due to its ability to form 
bound states (primarily with 4 He). When such bound 
states occur, they catalyze certain nuclear reactions 
such as 4 He(D, 7) 6 Li which is normally highly sup- 
pressed due to the production of a low energy 7 whereas 
the bound state reaction is not [55]. In Fig. [SJi, the 
( m i/2: m o) plane is displayed showing explicit element 
abundance contours [50] when the gravitino mass is 
m.3/2 = 0.2mo in the absence of stau bound state ef- 
fects. To the left of the solid black line the gravitino 
is the not the LSP. The diagonal red dotted line cor- 
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Fig. 8. Examples of mSUGRA (mi/s,mo) planes with 
contours o/tan/3 superposed, for /i > and (a) the simplest 
Polonyi model with Ac/mo = 3— \/3, and (b) Ao/mo = 2.0, 
all with Bo — Ao — mo. In each panel, we show the re- 
gions excluded by the LEP lower limits on MSSM particles 
and those ruled out by b —* S7 decay (medium green shad- 
ing): the regions favoured by g^ — 2 are very light (yellow) 
shaded, bordered by a thin (black) line. The dark (chocolate) 
solid lines separate the gravitino LSP regions. The regions 
favoured by WMAP in the neutralino LSP case have light 
(turquoise) shading. The dashed (pink) line corresponds to 
the maximum relic density for the gravitino LSP, and re- 
gions allowed by BBN constraint neglecting the effects of 
bound states on NSP decay are light (yellow) shaded. 



responds to the boundary between a neutralino and 
stau NSP. Above the line, the neutralino is the NSP, 
and below it, the NSP is the stau. Very close to this 
boundary, there is a diagonal brown solid line. Above 
this line, the relic density of gravitinos from NSP decay 
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is too high, i.e., 

^^-Q NS ph 2 > 0.12. (12) 

m NS p 

Thus we should restrict our attention to the area below 
this line. 

The very thick green line labelled 7 Li = 4.3 cor- 
responds to the contour where 7 Li/H = 4.3 x 10~ 10 , 
a value very close to the standard BBN result for 
7 Li/H. It forms a 'V shape, whose right edge runs 
along the neutralino-stau NSP border. Below the V, 
the abundance of 7 Li is smaller than the standard 
BBN result. However, for relatively small values of 
m i/2i the 7 Li abundance does not differ very much 
from this standard BBN result: it is only when ^ 
3000 GeV that 7 Li begins to drop significantly. The 
stau lifetime drops with increasing ttt, 1 / 2 7 and when 
r ~ 1000 s, at rrii/2 ~ 4000 GeV, the 7 Li abundance 
has been reduced to an observation- friendly value close 
to 2 x 10 -10 as claimed in [51] and shown by the (un- 
labeled) thin dashed (green) contours. 

The region where the 6 Li/ 7 Li ratio lies between 
0.01 and 0.15 forms a band which moves from lower left 
to upper right. As one can see in the orange shading, 
there is a large region where the lithium isotopic ratio 
can be made acceptable. However, if we restrict to D/H 
< 4.0 x 10 -5 , we see that this ratio is interesting only 
when 7 Li is at or slightly below the standard BBN 
result. 

Turning now to Fig. [9]b, we show the analogous 
results when the bound-state effects are included in 
the calculation. The abundance contours are identical 
to those in panel (a) above the diagonal dotted line, 
where the NSP is a neutralino and bound states do 
not form. We also note that the bound state effects on 
D and 3 He are quite minimal, so that these element 
abundances are very similar to those in Fig. [9^. How- 
ever, comparing panels (a) and (b), one sees dramatic 
bound-state effects on the lithium abundances. Every- 
where to the left of the solid blue line labeled 0.15 
is excluded. In the stau NSP region, this means that 
m i/2 ^ 1500 GeV. Moreover, in the stau region to the 
right of the 6 Li/ 7 Li = 0.15 contour, the 7 Li abundance 
drops below 9 x 10~ n (as shown by the thin green dot- 
ted curve). In this case, not only do the bound-state 
effects increase the 6 Li abundance when mj/2 is small 
(i.e., at relatively long stau lifetimes), but they also de- 
crease the 7 Li abundance when the lifetime of the stau 
is about 1500 s. Thus, at (m 1/2 ,m ) ~ (3200,400), 
we find that 6 Li/ 7 Li ~ 0.04, 7 Li/H ~ 1.2 x 10~ 10 , 
and D/H ~ 3.8 x 10 5 . Indeed, when m 1 / 2 is between 
3000-4000 GeV, the bound state effects cut the 7 Li 
abundance roughly in half. In the darker (pink) re- 
gion, the lithium abundances match the observational 
plateau values, with the properties 6 Li/ 7 Li > 0.01 and 
0.9 x 10~ 10 < 7 Li/H < 2.0 x 10~ 10 . 

Next we return to an example of an mSUGRA 
model based on the Polonyi model for which Aq /mo = 
3 — \/3 with the condition that m 3 / 2 = m . In Fig. [TUk . 
we show the mSUGRA model without the bound states. 
In the upper part of the plane, we do not have gravitino 
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Fig. 9. An(m 1 /2,nio) planes for Ao = 0, fi > and 
tan/3 = 10 with 7713/2 = 0.2mo (a) without and (b) with the 
effects of boundstates included. The regions to the left of the 
solid black lines are not considered, since there the gravitino 
is not the LSP. In the orange (light) shaded regions, the dif- 
ferences between the calculated and observed light- element 
abundances are no greater than in standard BBN without 
late particle decays. In the pink (dark) shaded region in 
panel b, the abundances lie within the ranges favoured by 
observation. The significances of the other lines and con- 
tours are explained in the text. 

dark matter. We see that 3 He/D 52 eliminates all but 
a triangular area which extends up to mo = 1000 GeV, 
when m 1/2 = 5000 GeV. Below the 3 He/D = 1 con- 
tour, D and 7 Li are close to their standard BBN val- 
ues, and there is a substantial orange shaded region. 
We note that 6 Li is interestingly high, between 0.01 
and 0.15 in much of this region. 

As seen in Fig. [T0b . when bound-state effects are 
included in this mSUGRA model, both lithium isotope 
abundances are too large except in the extreme lower 
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right corner, where there is a small region shaded or- 
ange. However, there is no region where the lithium 
abundances fall within the favoured plateau ranges. 
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Fig. 10. As in Fig. [H for the mSUGRA case with 
Ao/mo = 3 - V3. 

The prospects for direct detection in the mSUGRA 
models (when the neutralino is the LSP) is very similar 
to that of the CMSSM but slightly more difficult. In 
a plot such as that shown in Fig. the mSUGRA 
predictions are down by a factor of about 2 or less. 



9 The NUHM 

In the NUHM, the Higgs soft masses are treated inde- 
pendently from mo. In effect, this allows one to choose 
H and the Higgs pseudoscalar mass, ttta freely (up to 
phcnomenological constraints) . Two examples of planes 
in the NUHM are shown in Fig. [TTJ In panel (a), an 



m i/2~ m o plane is shown with ii = 700 GeV and ttia = 
400 GeV fixed across the plane [T8J . As usual, the light 
(turquoise) shaded area is the cosmologically preferred 
region. There is a bulk region satisfying this preference 
at m 1/2 - 150 GeV to 350 GeV and to ~ 100 GeV. 
The dark (red) shaded regions are excluded because a 
charged sparticle is lighter than the neutralino. As in 
the CMSSM, there are light (turquoise) shaded strips 
close to these forbidden regions where coannihilation 
suppresses the relic density sufficiently to be cosmolog- 
ically acceptable. Further away from these regions, the 
relic density is generally too high. At small TO1/2 and 
too the left handed sleptons, and also the sneutrinos, 
become lighter than the neutralino. The darker (dark 
blue) shaded area is where a sneutrino is the LSP. 

The near-vertical dark (black) dashed and light 
(red) dot-dashed lines in Fig. [11] are the LEP exclu- 
sion contours m x ± > 104 GeV and rrih > 114 GeV 
respectively. As in the CMSSM case, they exclude low 
values of mi/2, and hence rule out rapid relic annihila- 
tion via direct-channel h and Z° poles. The solid lines 
curved around small values of m 1 / 2 and too bound the 
light (pink) shaded region favoured by a M and recent 
analyses of the e + e~ data. 

A striking feature in Fig. [TlT a) when mi/ 2 ~ 450 
GeV is a strip with low fi x h 2 , which has bands with 
acceptable relic density on either side. The \cw-Q x h 2 
strip is due to rapid annihilation via the direct-channel 
A, H poles which occur when m x = ttia/2 = 200 GeV, 
indicated by the near-vertical solid (blue) line. These 
were found in the CMSSM [T0l[T2] . but at larger tan/3. 
At higher m 1 / 2 ~ 1300 — 1400 GeV, there is another 
region of acceptable relic density. This band, the tran- 
sition band, is broadened because the neutralino ac- 
quires significant Higgsino content as m^ 2 becomes 
greater than ii, and the relic density is suppressed by 
the increased W + W~ production. To the right of this 
band, the relic density falls below the WMAP value. 

Another example of an NUHM plane is shown in 
panel (b) of Fig. [TTl representing an ii — ttta plane 
with a fixed choice of values of tan [3 = 35, m 1 / 2 = 500 
GeV, m = 1000 GeV and A = [38]. Line types 
and shading are as in panel (a). The region with ac- 
ceptable relic density now takes the form of a narrow 
strip of values of /i ~ ±300 — > 350 GeV where the relic 
density lies within the WMAP range for almost all val- 
ues of rriA- The exception is a narrow strip centred on 
Tji a ~ 430 GeV, namely the rapid-annihilation funnel 
where m x ~ which would be acceptable if there 

is some other source of cold dark matter. Note the in- 
creased importance of the b — > sj constraint which 
now excludes almost all of [i < 0. 

By using the properties found in the ix — ttla plane 
discussed above, one can construct a ttia — tan (3 plane 
such that (nearly) the entire plane respects the cosmo- 
logical relic density constraint [3511^3] . These plane can 
be used as a particularly useful benchmark for existing 
and future searchs for Higgs bosons |53j . An example of 
such a plane is shown in Fig. 1121 with fixed TO1/2 = 500 
GeV and too = 1000 GeV. As one can see, the light 
shaded (turquoise) region now fills up the plane ex- 
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Fig. 11. Projections of the NUHM model on (a) the 
(mi/2, mo) planes for tan/3 = 10 and /i = 700 Gel/ and 
mA = 400 GeV and (b) the \i — mA plane with vti\n = 500 
GeV, m = 1000 GeV, tan/3 = 35 and A = 0. The (red) 
dot-dashed lines are the contours mh = 114 GeV as cal- 
culated using FeynHiggs [2^1, and the near-vertical (black) 
dashed lines are the contours m x ± = 104 GeV. The light 
(turquoise) shaded areas are the cosmologically preferred 
regions. The dark (brick red) shaded regions is excluded be- 
cause a charged particle is lighter than the neutralino, and 
the darker (dark blue) shaded regions is excluded because 
the LSP is a sneutrino. In panel (a) there is a very small 
medium (green) shaded region excluded by b — > s'y, at small 
mi/2- In panel (b), it excluded most of /j, < 0. The regions 
allowed by the E821 measurement of at the 2-a level, 
are shaded (pink) and bounded by solid black lines. 



BR(B U — > tv t ) and AMb s was performed recently in 
[55] . The lowest \ 2 value in this plane is Xmin = ^-4 
and the corresponding best-fit point is shown by a 
(red) cross as are the Ax 2 = 2.30 and 4.61 contours 
around the best-fit point. The (black shaded) region 
in the lower left is excluded at the 95 % C.L. by the 
LEP Higgs searches in the channel e + e~ — > Z* — > 

zh,H m\. 



m 1/2 = 500 , m = 1000 
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Fig. 12. The (m A , tan/3) plane for the NUHM bench- 
mark surface displaying the contours of A\ 2 found in a 
recent global fit to EWPO and BPO \35f . Also displayed 
is the expected sensitivity of present direct searches for the 
scattering of dark matter particles. 



Also shown in Fig. [TTj is the current reach of the 
XENON10 experiment [40]. The solid line in the figure 
correspond to the XENON10 bound obtained assum- 
ing U^n — 45 MeV. The corresponding reach obtained 
if £„n = 64 MeV is shown by the dotted curve. As one 
can see, in this case, the best fit point would already 
have been probed by the XENON10 experiment. If the 
strangeness content of the proton were zero, it would 
correspond to S^ N = 36 MeV, and the XENON10 
reach in this case is displayed by the dashed curve. 
Finally, we note that a future sensitivity to a cross 
section of 10 -8 pb would cover the entire surface. 

As can be inferred from the benchmark surface 
shown above, in the NUHM, current constraints al- 
ready exclude many interesting models [38] . These mod- 
els generically have large scattering cross sections. Scat- 
ter plots of showing the reach of CDMS and XENON10 
for both ZVjv = 45 MeV and 64 MeV are shown in Fig. 
1131 In comparison with Fig. [5] one sees that many more 
viable models are already currently being probed. 



cept for a vertical swath where s-channel annihilation 
through mA drives the relic density to small values. 
A likelihood analysis of this NUHM benchmark sur- 
face, including the EWPO M w , sin 2 9 eS , P Zl (g - 2) M 
and m h and the BPO BR(6 ->■ si), BR(B S -> fx + fj,~), 



10 A Hint of Higgs 

The CDF and DO Collaborations have already estab- 
lished important limits on the (heavier) MSSM Higgs 
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Fig. 13. Scatter plots of the spin-independent elastic- 
scattering cross section predicted in the CMSSM (a) for 
SirN = 45 MeV (b) and 64 MeV. 



bosons, particularly at large tan/3 [Ml 5 5. 5 6, 5 TUB'S]. 
Recently the CDF Collaboration, investigating the chan- 
nel 



PP 



h, H, A) 



(13) 



has been unable to improve these limits to the extent of 
the sensitivity expected with the analyzed integrated 
luminosity of ~ 1 fb _1 [57], whereas there is no in- 
dication of any similar effect in DO data [58] .Within 
the MSSM the channel (|13| is enhanced as compared 
to the corresponding SM process by roughly a fac- 
tor of tan 2 (3 /((l + A b ) 2 + 9) [59], where A b includes 
loop corrections to the 4>bb vertex (see [59] for details) 
and is subdominant for the t + t~ final state. Corre- 
spondingly, the unexpected weakness of the CDF ex- 
clusion might be explicable within the MSSM if rriA ~ 
160 GeV and tan/3 > 45. 



In the CMSSM, there are no parameter sets con- 
sistent with all the experimental and theoretical con- 
straints. In particular, low values of are required 
to obtain low values of uia- For rriA — 160 GeV and 
the other parameters tuned to fulfill the B physics con- 
straints, there are no CMSSM solutions for which the 
lightest Higgs mass is not in direct contradiction with 
experimental limits [53]. However, we do find that all 
the constraints may be satisfied in the NUHM [SOllBT] . 

We consider first the (mi/2,mo) plane shown in 
panel (a) of Fig. [H which has tan/3 = 45, /u = 370 GeV 
and A = -1800 GeV, as well as M A = 160 GeV. As 
before, the dark (brown) shaded region at low mo is 
forbidden, because there the LSP would be the lighter 
stau. The WMAP cold dark matter constraint is satis- 
fied only within the lighter (turquoise) shaded region. 
To the left of this region, the relic density is too small, 
due to s-channel annihilation through the Higgs pseu- 
doscalar A. As m 1 / 2 increases away from the pole, the 
relic density increases toward the WMAP range. How- 
ever, as rai/2 is increased, the neutralino acquires a 
larger Higgsino component and annihilations to pairs 
of W and Z bosons become enhanced. To the right 
of this transition region, the relic density again lies 
below the WMAP preferred value. The shaded region 
here is therefore an overlap of the funnel and transition 
regions discussed in [T5]. The BR(_B S — > /i + /i~) con- 
straint [62] is satisfied between the outer black dash- 
dotted lines, labelled 10 -7 , representing the current 
limit on that branching ratio. Also shown are the con- 
tours where the branching ratio is 2 x 10~ 8 , close to 
the sensitivity likely to be attainable soon by CDF 
and DO. Between these two contours, there is a strong 
cancellation between the flavor- violating contributions 
arising from the Higgs and chargino couplings at the 
one-loop level and the Wilson coefficient counterterms 
contributing to BR(i3 s — > n + [63] . 

The dash-dotted (red) line shows the contour cor- 
responding to m/i = 114 GeV, and only the region 
to the right of this line is compatible with the con- 
straint imposed by m/j (it should be kept in mind 
that there is still a ~ 3 GeV uncertainty in the pre- 
diction of m/j [24]). Also shown in pink shading is the 
region favoured by (g — 2) p at the two-cr level. The 
one- and two-cr contours for {g — 2)^ are shown as el- 
liptical dashed and solid black contours, respectively. 
The region which is compatible with the WMAP relic 
density and mh, and is also within the two-cr (g — 2)u 
experimental bound, has BR(£? S — > > 2x 10 . 

The measured value of BR(& — > sj) is in agreement 
with the theory prediction only to the left of the solid 
(green) region. 

We see that there is a narrow wedge of allowed 
parameter space indHa), which has mi/2 ~ 600 GeV 
and m - 700 to 1100 GeV. The BR(& -> sj) con- 
straint is satisfied easily throughout this region, and 
(g — 2)^ cuts off the top of the wedge, which would oth- 
erwise have extended to mo S> 1500 GeV. Within the 
allowed wedge, mt is very close to the LEP lower limit, 
and BR(B S -> n + n~) > 2 x 10~ 8 . If M A were much 
smaller (< 130 GeV), there would be no wedge consis- 
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Fig. 14. The NUHM parameter space as a function 
of mi/2 and mo for /i = 370(400) GeV and Ao = 
-1800(-2000) GeV in panels a (b). We fix M A = 160 GeV, 
tan/3 = 45, m* = 171.4 GeV and m b = 4.25 GeV. For the 
description of the various lines and shaded areas, see the 
text. 



tent simultaneously with the i?cDM, witi and BR(_B S — > 
/i + /i~) constraints. 

Given the matrix element uncertainties for direct 
detection summarized above, we show in 114b ) the l-o~ 
lower limit on the calculated value of the elastic cross 
section as compared to the CDMS upper limit [3S] . In 
the portion of the plane to the left of the (orange) 
dotted line, the lower limit on the calculated spin- 
independent elastic cross section is smaller than the 
CDMS upper bound, assuming the canonical local den- 
sity. Whilst we have assumed E^n — 45 MeV, the cal- 
culated lower limit effectively assumes zero strangeness 
contribution to the proton mass, i.e., y = 0. In the 
region of interest, the lower limit on the calculated 
cross section is about 80% of the CDMS upper bound, 



whereas with a strangeness contribution of y — 0.2, 
the cross section would exceed the CDMS bound by a 
factor of ~ 3. Thus, if Nature has picked this corner of 
the NUHM parameter space, we expect direct detec- 
tion of dark matter to be imminent. Consistency with 
the XENON10 limit would further require a reduction 
in the local dark matter density (to its lower limit). 
Intriguingly, the XENON10 experiment has seen some 
potential signal events that are, however, interpreted 
as background. 

Panel (a) of [TBI explores the (/j,,Aq) plane for the 
choice (mi/a, mo) = (600, 800) GeV, values close to the 
lower tip of the allowed wedge in 114( a). In this case, 
the region allowed by the BK(B S — > /i + /i~) constraint 
is below the upper dash-dotted black line, and the LEP 
to/j constraint is satisfied only above the dash-dotted 
red line. We see that only a restricted range 360 GeV < 
\i < 390 GeV is compatible with the dark matter con- 
straint. This corresponds to the transition strip where 
the neutralino is the appropriate bino/Higgsino com- 
bination. To the left of this strip, the relic density is 
too small and to the right, it is too large. Only a very 
restricted range of Ao ~ —1600 GeV is compatible si- 
multaneously with the TOft and BR(£? S — ► /i + /i~) con- 
straints. Very large negative values of Ag are excluded 
as the LSP is the lighter stau. On the (/i, Aq) plane, the 
elastic scattering cross section is a rapidly decreasing 
function of \i and is almost independent of Aq. Indeed, 
NUHM points excluded by CDMS (or XENON10) gen- 
erally have low values of // and Ma E] ■ Values of /i > 
355 GeV are compatible with CDMS if the strangeness 
contribution to the proton mass is negligible. For this 
choice of parameters, the entire displayed plane is com- 
patible with BR(fo — * sj) and [g — 2) M . 

Panel (b) of Fig. ITSI shows what happens if TOi/ 2 
is increased to 700 GeV, keeping too and the other 
inputs the same. Compared to Fig. 115( a). we see that 
the WMAP strip becomes narrower and shifts to larger 
|U ~ 400 GeV, and that BR(6 — > 57) starts to exclude 
a region visible at smaller fj,. If TO1/2 were to be in- 
creased much further, the dark matter constraint and 
BR(fe — > s >) would no longer be compatible for this 
value of Too. We also see that, by comparison with 
Eta), the BR(_B S — * constraint has moved to 

lower Ao, but the mh constraint has dropped even 
further, and mh > 114 GeV over the entire visible 
plane. The net result is a region compatible with all 
the constraints that extends from Aq ~ —1850 GeV 

down to A Q 2150 GeV for /j, ~ 400 GeV, with 

a coannihilation filament extending to larger /1 when 

A 2200 GeV. Once again all of the WMAP strip 

in this panel is compatible with CDMS. 

The larger allowed area of parameter space is re- 
flected in panel (b) of Fig. El which has \i = 400 GeV 
and A = -2000 GeV, as well as tan (3 = 45 and M A = 
160 GeV as before. In this case, we see that a substan- 
tial region of the WMAP strip with to 1/2 ~ 700 GeV 
and a width #TOi/2 ~ 100 GeV, extending from too ~ 
750 GeV to higher too is allowed by all the other con- 
straints. The BR(_B S — > /i + ^i~) and mh constraints 
have now moved to relatively low values of TO1/2, but 
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Fig. 15. The NUHM parameter space as a function of \i 
and A for ra 1/2 = 600(700) GeV and m = 800 GeV in 
panels a (b). We fix M A = 160 GeV, tan/3 = 45, m t = 
171.4 GeV and nib = 4.25 GeV. For the description of the 
various lines and shaded areas, see the text. 



we still find BR(B S -> > 2 x 1(T 8 and m h 

close to the LEP lower limit. The (g — 2) M constraint 
truncates the allowed region at mo ~ 1050 GeV. Once 
again, the allowed region is compatible with CDMS 
(to the left of the (orange) dotted line) provided the 
strangeness contribution to the proton mass is small. 

Clearly our anticipation for the discovery of super- 
symmetry or new physics is great. If the CDF hint of 
the Higgs boson is realized, we should expect to see nu- 
merous departures from the standard model including 
the decay of B s to (J, + fi~, BR(b — ► sj) should show 
deviations from the standard model value, and dark 
matter should be discovered by CDMS or XENON10 
in the near future. Similarly if an observation of new 
physics in any one of these areas occurs, we should be 



ready for flood of new data showing departures from 
the standard model. 
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